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In one dimensional transport problems the scattering matrix S is decomposed into a block 
structure corresponding to reflection and transmission matrices at the two ends. For S a 
random unitary matrix, the singular value probability distribution function of these blocks 
is calculated. The same is done when S is constrained to be symmetric, or to be self dual 
quaternion real, or when S has real elements, or has real quaternion elements. Three meth- 
ods are used: metric forms; a variant of the Ingham-Seigel matrix integral; and a theorem 
specifying the Jacobi random matrix ensemble in terms of Wishart distributed matrices. 



1 Introduction 

In mesoscopic physics, the conductance of certain quasi one-dimensional wires containing scattering 
impurities exhibit the phenomenum known as universal conductance fluctuations. As first predicted 
theoretically ^ II lj . and soon after observed experimentally [TfJ, for such wires the variance of the 
conductance is of order (e 2 /h) 2 , independent of sample size or disorder strength. Furthermore the variance 
decreases by precisely a factor of two if time reversal symmetry is broken by a magnetic field. This 
conductance problem is fundamental for its relation to the Landauer scattering theory of electronic 
conduction (see 0] and references therein) , and to time reversal symmetry. It is further fundamental for 
its relation to random matrix theory. Let us revise these points by way of background and motivation 
for the specific problem of this paper. 

In the theoretical description of the conductance problem, basic quantities are the electron fluxes at 
the left and right hand edges of the wire. These are specified by an n-component vector I and an m- 
component vector I' specifying the complex amplitudes of the available plane wave states travelling into 
the left and right sides of the wire respectively, as well as an n-component vector O and an m-component 
vector O' for the same states travelling out of the left and right sides of the wire. For definiteness it will 
be assumed that n > m. 



By definition, the scattering matrix S relates the flux travelling into the conductor to that travelling 



out, 
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(1.1) 



The scattering matrix is further decomposed in terms of reflection and transmission matrices by 



S = 



(1.2) 
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Flux conservation requires 

\i\ 2 + \r\ 2 = \6\ 2 + \d'\ 2 , 

and this implies that S must be unitary. Furthermore, by relating S to an evolution operator and thus 
a Hamiltonian, one can argue (see e.g. [Hj) that S must be symmetric if the system has a time reversal 
symmetry with T 2 = 1 , and a self dual quanternion matrix when there is a time reversal symmetry with 
T 2 = -1. 

The immediate relevance of the above formalism is seen by invoking the Landauer scattering theory 
of electronic conduction. According to this formalism, the conductance G is given in terms of the 
transmission matrix t mxn (or t' nXm ) by the so called two probe Landauer formula 

G/G = Tt(th) = Tr{t'W) (1.3) 

where Go = 2e 2 /h is twice the fundamental quantum unit of conductance. Thus to compute G it suffices 
to know the distribution of the eigenvalues of tH (or t'H'). In fact the matrix S can be decomposed in a 
form which isolates these eigenvalues. 

For definiteness suppose there is no time reversal symmetry, so S is a general unitary matrix. Each 
block of S can then be decomposed according to a general singular value decomposition. For example, 
for the block r nxn we have 

r nxn = U r A r V r , 

where U r , V r are unitary matrices and A r is a rectangular diagonal matrix with entries equal to the positive 
square roots of the eigenvalues of r*r (these eigenvalues are between and 1 since r'r + t't = I nxn ). The 
unitarity of S inter-relates the matrices U r , U r >, . . . (see e.g. ^j]) an d implies the decomposition 
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(1.4) 



(1.5) 



Symmetries relate V r ,V r i to U r ,U r ,. As already remarked, in the case of a time reversal symmetry 
T 2 = 1, S must be symmetric, while for T 2 = — 1, S must be self dual quaternion real. These symmetries 
require that 

V} = Uj, V r \ = Ul and V} = U? ', V r \ = Ufi (1.6) 

respectively. In (|1.6|) the operation D, for annxn matrix A with real quaternion elements regarded as 
a 2n x 2n matrix with complex elements, is specified by 



A D = Z 2n A'Z^ 



T/7-1 



z 9 
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The relevance of random matrix theory comes from hypothesizing that apart from the constraint 
imposed by a time reversal symmetry, in the regime relevant to universal conductance fluctuations, the 
scattering matrix S is effectively a random matrix chosen with Haar (uniform) measure. With this 
assumption the non-zero elements of A t , which in turn are equal to the square root of the non-zero 
eigenvalues of tH and are thus in the interval (0,1), has a joint probability density function (p.d.f.) 
proportional to 



n A f n m-^. 



+ 1 - 1/j8, 



(1.7) 



l<j<k<m 
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where (3 = 1, 2,4 according to S being constrained to be symmetric, no constraints, or constrained to be 
self dual quanternion. In its full generality this result is stated without proof in Beenakker eq. (2.9)], 
where it is attributed to Brouwer. In the case n = m it was derived in With knowledge of (|1.7|l . 

it can be deduced from ((1.3|) that with a fixed |3] 

lim Var(G/G ) = (1.8) 

which is in quantitative agreement with experiment. 

However the issue of fluctuation formulas is not our concern here. Rather we seek to specify three 
distinct derivations of (|1.7|) in the case (3 = 2, each of which we believe to have particular features of 
interest. The first derivation to be given uses the method of metric forms. This method applies to all 
three cases, and in fact details will be given in the case (3 = 1. The second derivation uses versions of 
the Ingham-Seigel integral; while the third relies on knowledge of the Jacobi random matrix ensemble 
as derived from Wishart distributed matrices. These latter two derivations generalize (|1.7(l in the case 
(3 = 2, and also allow (3 generalizations relating to the block decomposition of unitary matrices, however 
now with (3=1 refering to the elements of S being real, and (3 = 4 to the elements being real quaternion. 
Some aspects of statistical properties of the block decomposition of real orthogonal matrices (i.e. unitary 
matrices with real elements) have been previously given in |18| . 



2 Metric forms 



Let X = [xjk]j.k=i,...,N be an N x N matrix. Let {x^}, where fi labels both rows and columns, be the 
set of independent real and imaginary parts in X. The metric form of the line element ds is defined by 

(ds) 2 := Tr(dXdX^) = ( dx n) 2 + 2 H (dx^) 2 . 

}jl diag. 



M upper 
triangula 



The corresponding volume measure is 



(dX)= /\ dx^ f\ dx». 

fj, diag. 



(2.1) 



i± upper 



Suppose now a change of variables {x^} i— > {y^} is made such that (ds) 2 is a symmetric quadratic form 
in {dy^}, 

(ds) 2 = Y 9^^dy^dy„, 5^ = g v ^. 



The corresponding volume measure is then (see e.g. |12j) 

(det^D 1 / 2 /\d Vli 



(2.2) 



In [Hj this formalism has been used to derive the result 1)1.7(1 in the case (3 = 2. Here the details will 
be given for the case (3 = 1, when 11.4(1 reads 



(2.3) 
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will be given. 



3 



For a general matrix A, let SA := A'dA, where dA denotes the matrix of differentials of the elements 
of A. Using the fact that for A unitary (SA)^ = —SA, it follows from l|2.3[l that 



(2.4) 



In Ti(dSdS^), the right hand side of (|2.4() can effectively be substituted for dS. Doing this and simplifying 
using 

r 5U r 

5U r > 
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5U r > 







SU r < 



Tr 



LUL = 



gives 



TrfUSU^) = Tr AiAi + Tr(A 2 + A 2 ) + Tr(A 3 A 3 ) + Tr(dLdtf) 
where, with / denoting the identity matrix, 



(2.5) 



Ax 
A 2 
A 3 



I - AtAfSUr - 5Ur\JI - A t Aj 

{A t 5U r > - SU r At)(AjSU r - SUr'Aj) 



I - A t Aj 5U r , - 5U r >\/l - A t AT 



We recall that the diagonal elements of A t are the positive square roots of the eigenvalues of tH. 
Because t = t mxn and thus has rank m, t't must have n — m zero eigenvalues, so with the diagonal 



elements of A t Aj denoted Xj (j = 1 



7 " ' " 5/5 



as is consistent with the notation of Q1.7[l . we have X m +i 



■ ■ ■ = X n — 0. Using this fact we see that in component form the four terms in (|2.5|l . T\, . . . , T4 say, can 
be expanded to read 



Ti = J2(l-Xl)\(SU r ) kk -(SU r ) kk \ 2 + \{^^+^ 



fe=i 



l<fc<2<n 



A? ) \(SU r ) lk - (SU r )l k f 



m n 



-( E +E E )\ 

l<k<l<m k=ll=m+l 



1-Xf-Jl-X 2 k ) \{8U r )i k + {8U r ) lk \< 



To 



J2 X l\( dU r')kk~(SU r ) kk \ 2 + {^( X l+ X k) 2 \(SU r ') lk ~{SU r ) lk \ 
k=l l<k<l<m 

tci n 



+-(x l -x k ) 2 \(su 1 ,) lk + (su r ) lk \ 2 } + J2 E x l\(su r : 



Ik 



k=l 2=m+l 



T 3 = ^(l-A2)|( ( 5C7 I .0 fefe -(<5{7,0fefe| 2 + E 2(V 1_A " + V 1_A *) \i 5tj r')ik-{SU r ')ik\ 



k=l 
1 



l<fc<Z<m 



T 4 



i=i 1 



A r 



In general a symmetric unitary matrix has the same number of independent elements as a real sym- 
metric matrix of the same rank, so S has \{ra + n)(m + n + 1) independent elements. We thus seek this 
same number of independent differentials in Ti, . . . , T4. These are 



(SU r ){% \<k<n («7 r 0j& 1 < k < m 



(i) 
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and 

as well as 
and 



(SU r )\l\ (5U r )\l (SU r )\l\ (5U r )\l 1 < k < I < m, 
(5U r )$ , (6U r )$ , l<k<mkm + l<l<n 



{SU r ) 



(i) 
Ik • 



1<1 <k < 



d\j 1 < j ; < m, 



which indeed tally to \(rn + n){m + n + 1). 

We see from the expressions for T 1; T 2 , T 3 that the contribution to the metric form from the differentials 
with subscripts Ik such that 1 < k < I < m is 



l<k<l<m 



(i)\2 
Ik 



i-Ji-xfJi-xt 



] )1 + (5UP)1 



-2XiX k {5U^) lk {5U^) lk + 2X l X k (SU i ; ) ) lk (SUj:% k . 



(2.6) 



Setting 1 — A 2 Wl — A 2 , =: a for notational convenience, this portion of the metric form contributes to 
(det^fc]) 1 / 2 in l|2.2[l 2x2 block factors which is proportional to 



n 

k<l 



l + a -A;Afc 
-A;A fc l + a 



1/2 



1 - a A;A fc 
AjAfe 1 - a 



1/2 



(2.7) 



For m + l<fe<i<n, the coefficient of (#^r)ifc in T\ is independent of the A's and so for the present 
purposes can be ignored. For 1 < k < m, m + 1 < I < n, we see from Ti and T 2 that the corresponding 
contribution to the metric form is 



m n 



4 E E (i + v/i - ad((^)^) 2 + (i - V 1 - ^)((jc/ r )^) 2 

k=l l=m+l 

This contributes to the volume form a factor proportional to 



m n 



n n **=n* 

k=l l=m+l 



n—m 
k 



(2.8) 



k=l 



We read off that the contribution to (detff^-fe]) 1 ' 2 from the coefficients of the terms (dAj) 2 in T4 is 



proportional to 



(2.9) 



It remains to calculate the contribution from the differentials on the diagonal. For m + 1 < k < n, 
the coefficient of (SU r ) k % k is a constant so these differentials can be ignored. We read off from Ti,T^,T^ 
that the contribution to the metric form from the remaining differentials is 

ra 

2^(2 - Xt)(((SU r ) k ir + (6U r ,)i) 2 ) + 2Xt(SU r ) k V l(SU r ,) k l 
fe=i 

The contribution to the volume form is thus proportional to 

1/2 



n 



A 2 



XI 



oc 



nc 1 -^) 1 



/2 



(2.10) 



k=l 



Multiplying together l|2.7jl - H2.10j) gives (|1.7jl in the case f3 = 1. 
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3 Matrix integrals 



Let 



4 2 l(Q™) := / e^ H ^ (det(iT m - nl m )) % \dH m ) (3.1) 

where H m ,Q m are to x to Hermitian matrices, and suppose n > to, Im/i > 0. For Q OT positive definite, 
it has been proved by Fyodorov 7 that 

m(m+l)/2-m(_i\m(m-l)/2 . ■ s n-m . 

= ^-i J \ (d<Q m ) e i^«-. (3.2) 



nr="or(n-i) V ^2 

(Here we have taken (dH m ) := Iljli Y\j<k ^jk^jk wmcn differs by a factor of 2 in the product 
over j < fc to the convention adopted in [7].) This matrix integral may be regarded as being of the type 
first evaluated by Ingham and Siegel (see [7] and references therein). In a subsequent work 1)3. lfl was 
used to derive l|1.7l) in the case (3 — 2. Here we will show that this derivation can be used to derive a 
generalization of 1|1.7[) in the case (3 = 2, and this generalization can be further extended to the cases 
0=1 and = 4, using suitable variants of (|3.1|l . where now (3 = 1 and 4 refers to the decomposition 
<|1.2H with S having real and real quaternion elements respectively. 

Define by Im}n(,Q m ) the matrix integral (|3.1|l with H rn and Q m now real symmetric. Also, define by 
Im}n(Qm) the same matrix integral but with H rn and Q m now self dual quaterion Hermitian matrices 
(such matrices regarded as 2m x 2to Hermitian matrices are doubly degenerate; adopt the convention 
that the operations Tr and det include only distinct eigenvalues). In [7] it is remarked that the method 
of derivation given therein to deduce that the evaluation of Im'n(Q m ) can also be used to deduce the 
evaluation of I^niQm), which reads 



m,n/2 ' 



nm m(m+3)/2 -m(m+l)/2 . a , (n— m-l)/2 

(Q m ) = — 1 fdet(-Q m )V e *^0». (3.3) 

nr=o r((n-i)/2) V { 2 Wm) J y ' 

Applying the same method to 4;„(Q m ) gives 

... (9iY n Tr m2 / 7 \ 2(n— m+1/2) 

I ( X(Qm) = * .M ^Qm)) } e*^«-. (3.4) 

LIj=o r ( 2 (n- j)) v 2 / 

Hence, all three cases we have 

j£i»/2«m) = (det(Q m ))^)(-™ +i-V« 6 ^«3. > (3.5) 

where Cm,„ is independent of Q m . 

Let U be an N x N random unitary matrix, with real {[3 = 1), complex {(3 = 2) and real quaternion 
(,9 = 4), chosen with Haar measure. Generalizing 1)1. 2[1 . decompose J7 into blocks 



C7 = 

Since J7 is unitary, we require 



^-nixri2 Cnix(A r — n 2 ) 

B(N-m)xn 2 -D(AT-ri 1 )x(Ar-ri 2 ) 



(3.6) 



i4At + CCt = J B1 , (3.7) 

together with three similar equations involving B and £> which given A and C can always be satisfied. 
The idea of ^S] is to regard (|3.7I) as a constraint in the space of general rectangular matrices A, C with 
entries of the type required by the index (3. Thus in this viewpoint the distribution of A is given by 

6(AA* + CC^ - I n . 2 ){dC). (3.8) 



G 



In (|3.8fl the delta function is a product of scalar delta functions, one for each independent real and 
imaginary component of A' A + C^C — I n2 . It is proportional to the matrix integral 

e -iTr(H(AAl+CCl-I„ 2 ))( dH) (39) 

where H is an Hermitian matrix with elements of the type f3. Following j8| we would like to substitute 
<|3.9|l for the delta function in Q3.8[l . and change the order of integration. The integrations over (dC) are 
simply Gaussian integrals. For the resulting function of H to be integrable around H = 0, the replacement 
H H — i[il ni in the exponent of l|3.9|) must made. Doing this and computing the Gaussian integrals 
gives that (|3.8|) is proportional to 

lim J ( det (H iM 7 ni )) "'^V^W* -a* )) {dH) 

Evaluating the matrix integral using (|3.5|) shows that the distribution of A is proportional to 

^det(7 ni - AA^)J . (3.10) 

Note that for this to be normalizable, we must have 

N-m-n 2 >0. (3.11) 

Suppose that in addition to l|3.11(l we have n\ > ri2- Then A has rank n\ and so AA^ has n\ — ri2 
zero eigenvalues and l|3.10|l can be written 

( det(J„ 2 -A^A)\ . (3.12) 
Setting Y ~ A 1 A we know (see e.g. 6 ) that 

dA oc (detYyP/^-^-WidY), (3.13) 
so we have from (|3.12f) that the distribution of Y is proportional to 

(detF)(' 3 / 2 )^^+ 1 - 2 /«)(det(/„ 2 -y)) (3.14) 



Denote the eigenvalues of Y by y\, . . . , y ri2 . Using the fact that the eigenvalue dependent portion of the 
Jacobian for an Hermitian matrix, with elements of the type f3, when changing variables to its eigenvalues 
and eigenvectors is Ylj <k \yk — 2/j l' 3 , we read off from l|3.14|l that the eigenvalue p.d.f. of Y is proportional 
to 

ri2 "2 

JJ I/ (/J/2)(nx-« 3+1 -a//J)( 1 _ 1/i) aJ/2)(JV-nx-n 3 +l-a//J)) FJ \y k - Vj f . (3.15) 

j=l j<k 

In the case 

n\ = 712 = m, N — n + m, (3.16) 

the matrix A in (|3.t)|l coincides with r in l|1.2|l . and the yj in l|3.15|l are then the eigenvalues of r'r. The 
non-zero values of the singular values of the submatrix t in (|1.2|l are the Aj's in (|1.7|l . The decomposition 
(|1.4|l tells us that yj = 1 — Xj. Indeed, making this change of variable in (|3.15|l . and making the 
substitutions lf>|) . reclaims (|1.7(l in the case /3 = 2. 
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For (3 = 2 and general ri\ > n 2 , N > n\ + n 2 , a result equivalent to (|3.15[) is derived in a recent work 
of Simon and Moustakas JS]. Motivated by a quantum dot problem with three leads, they decomposed 
the N x N scattering matrix S into a 3 x 3 block structure 





ru 


tl2 


*13 


3 = 


tix 


'"22 






hi 


^32 


h:i 



where each r« is Ni x Ni and Uj is Ni x Nj, with TV = iVi + N2 + N3 . For S a random unitary matrix with 
Haar measure, it is shown that the eigenvalues of t\ 2 ti 2 have p.d.f. given by (|3.15|l with (3 = 2, ni = N, 
TI2 = N 2 . This is consistent with our result because the distribution of S is unchanged by interchanging 
rows and columns. The first and second block columns in l|3.17|l can be interchanged, effectively giving 
the decomposition i|3.6[) with i 12 = An iX n 2 . 

We remark that for yj <C 1 the term involving (1 — j/j) in l|3.15|l can be ignored. The resulting p.d.f. is 
then the Xj <C 1 limit of the singular values p.d.f. of the matrix product X'X, where X is an n\ x n 2 
random Gaussian matrix with real (/? — 1), complex (/3 — 2) or real quanternion elements (/? = 4) (see 
e.g. [S]). This is consistent with a recent result of Jiang ^U], who quantifies the degree to which the 
entries of orthogonal or unitary random matrices can be approximated by entries of Gaussian random 
matrices with real or complex entries respectively. 



4 Random projections 

In a recent work Collins computed the distribution of certain products of random orthogonal projec- 
tions. Here this approach, suitably modified, will be used to derive (|3.15l) . Knowledge is required of the 
following result |131 |S] 

Proposition 4.1. Let c and d be n\ x m and n 2 x m Gaussian random matrices, all elements identically 
and independently distributed (i.i.d.) where n 2 ,n 2 > m, indexed by the parameter (3 = 1,2 or 4. The 
parameter (3 specifies that the elements are real {(3 = 1), complex {(3 = 2) or real quaternion (j3 = 4). Let 
C = c+c, D = dU. The distribution of J := (C + D)- l l 2 C(C + D)~ 1 / 2 is proportional to 

(det J) Q/3 / 2 (det(/ - J)) b ^ 2 (4.1) 

with 

a = ni-m+l-2//3, b = n 2 - m + 1 - 2/j3. 

The p.d.f. H4.1|) is referred to as the Jacobi ensemble (see e.g. To make use of Proposition 14. II 
let X be a TV x to, N > m i.i.d. Gaussian random matrix indexed by [3. Define Q\ to be the top 
n x N submatrix of an N x N unitary matrix indexed by (3. Because QiQ\ = I n , we have that Q\X is 
distributed as an n x to i.i.d. Gaussian random matrix indexed by f3 (see Thm. 2.3.1]). Let Q 2 be the 
bottom (N — n) X N submatrix of U so that 

Q\Qi + QlQ2 = In- (4.2) 

Because Q2Q2 — ^N-n, Q2X is distributed as an (N — n) x to i.i.d. Gaussian matrix indexed by (3. 
Now set 

Y = X^Q\Q X X, Z = X^Q\Q 2 X. (4.3) 
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Denoting by Wp^\n) the (Wishart) distribution of the matrix product A^A, where A is an n xp, (n > p) 
i.i.d. Gaussian random matrix indexed by (3, we see that the distribution of Y is Wm\n) while the 
distribution of Z is Wm\N — n). The result of Proposition ^. II tells us that the eigenvalue p.d.f. of 

(Y + Z)-V*Y(Y + Z)- 1 ' 2 = (X^Xy^YiX^X)- 1 / 2 (4.4) 

is given by the Jacobi ensemble (|4.1() with n\ = n, n 2 = N — n. 
According to the singular value decomposition, we can write 

X = E/1AU2 (4.5) 

for Ui an N x N unitary matrix, U 2 an m x m unitary matrix, and A an N x m diagonal matrix, with 
diagonal entries equal to the positive square root of the eigenvalues of X^X. It follows from (|4.5|l that 

{X^X) 1 ' 2 = U 2 [K\ m , m U 2 (4.6) 

where [A] mjjn refers to the top in x m sub-block of A, while l|4.6|l and l|4.5|l together imply 

X^X)- 1 ' 2 = [U x ] N , m ul (4.7) 

Substituting for Y in (|4.4|) according to (|4.3fl . then using (|4.7|l to substitute for X(X^X)- 1 / 2 tells us 
that the distribution of 

U2[Ux] ] N , m Q\Qx[U x ] N , m Ul (4.8) 

is given by (|4.1|) . We note 

f/ 1 t [C/i]jv,m = diag(l,...,l,0,...,0) 

where there are m l's. Since the distribution of Q\ is unchanged upon multiplication by a unitary matrix, 
the distribution of 1)4. 8fl is the same as the distribution of 

U 2 [Q{Qi} m , m ul (4.9) 
In terms of the decomposition (|3.6f) this in turn has the same distribution as 



A*A ni _ . (4.10) 

Hence we have the distribution of 14.10J1 is given by 14.1J1 with 

a = n - m + 1 - 2//3, b = N - n - m + 1 - 2/0, 
which is in precise agreement with (|3.14[) upon the replacements ri\ 1— > n, n 2 ^ m in the latter. 
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